Origin of pressure induced second superconducting dome in A y Fe2- x Se 2 [^4=K, (Tl,Rb)] 
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Recent observation of pressure induced second superconducting phase in A v Fq2- x ^2 [A=K, (Tl,Rb)] calls 
for the models of superconductivity that are rich enough to allow for multiple superconducting phases. We 
propose the model where pressure induced changes in Fermi surface topology lead to appearance of the hole 
pocket at F point at high pressure. We develop low-energy effective model, derived from first-principles band- 
structure calculation, to suggest the phase assignment where low pressure superconducting state with no hole 
pocket at F point is a nodeless d-wave state. It evolves into a s ± state at higher pressure when the Fermi 
surface topology changes and hole pocket appears. We analyze the pairing interactions using five band tight 
binding fitted band structure and find strong pairing strength dependents on pressure. We also evaluate the 
energy and momentum dependence of neutron spin resonances in each of the phases as verifiable predictions of 
our proposal. 

PACS numbers: 74.70.Xa, 74.20.Rp, 74.25.Jb, 74.40.-n 



High— T c superconductivity often occurs when the system 
is driven from its pristine phase to the verge of magnetic quan- 
tum critical point via external controlled parameters such as 
chemical doping, magnetic field or pressure (P) in most of 
the cuprates, heavy fermions, pnictides, and organic super- 
conductors. However, several recent breakthrough discover- 
ies of a second superconducting (SC) dome — completely iso- 
lated or slightly connected to the first SC dome — without 
the intervention of any competing order as a function of P 
in several families of high-T c superconductors, GH3 and/ or 
extreme chemical doping in KFe2As2,|4] LaFeAsOi _9H^t5l 
or strainO have questioned this expectation. Much higher 
optimal T c value than that of the first dome, as well as a sub- 
stantial increase in electronic mass with P as obtained in iron- 
based compounds suggest an interesting and exotic phenom- 
ena of superconductivity along this tuning axis. It is notable 
that in stochiometric SrFe2As2, a SC dome appears both as a 
function of P,|7 ] as well as crystallographic strain (6l near the 
quantum critical point of spin-density wave, as often observed 
in other high-T c superconductors, further supporting the no- 
tion that the P induced superconductivity is unconventional. 
In this letter, we present a model that allows us to capture the 
onset of the second SC phase with different pairing symmetry 
as a function of P. 

Our approach is based on the weak or intermediate cou- 
pling scenario in which the shape of the Fermi surface (FS) 
topology plays a key role in creating pairing instability at the 
'hot-spot'.[8| For such case Cooper pairing arises from re- 
pulsive interaction with sign-reversal pairing symmetry con- 
strained by the FS topology and crystal symmetries. This the- 
oretical framework consistently describes d— wave pairing in 
cuprates. 1H1 and Ce-based heavy fermions,[9] s ± -pairing in 
iron-pnictides and chalcogenides,| 10 ] and nodeless d-wave in 
A y Fe 2 -xSe 2 [A=K,Cs, Rb,(Tl,Rb), (T1,K)] families. EHH 

Here, we focus on the latter family in which a second 
isolated dome is observed, as shown Fig. [TJa) (reproduced 
from Ref. 14). In an analogous 122 compound BaFe2As2, 




FIG. 1. (Color online) (a) Experimental phase diagram of T c vs. P 
for several A v Fq2- x Sq2 systems, reproduced from Ref. [3] (b) Mass 
renormalization for hole-pocket (a-FS) and electron pocket (/3-FS) 
for a related iron-pnictide compound BaFe2As2, deduced from quan- 
tum oscillation measurements (symbols). 1 14 1 The experimental data 
are normalized to their values at P = 0. Theoretical results of mass 
renormaliztion factor Z = m/m* for A y Fz2- x Sz2 compound (see 
text) plotted in solid line, (c) Ab-inito band structure of KFe2Se2 
at two representative Ps using TB-LMTO method (see text). Inset: 
Bands at P = (solid line) is multiplied by Z — 0.9 which matches 
well with that for P=12 GPa, affirming the fact that P increases ef- 
fective mass in this system as seen in experiment in (b). 



it is established that the effective band mass, m*, gradu- 
ally increases as a function of P, see the plot of 1/ra* in 
Fig. [TJb). Constrained by these experimental facts we pos- 
tulate that uniform P renormalizes the electronic structure, 
as also shown in first-principles calculations. fT5l and thereby 
affects the FS topology. In A v ¥q2-x^2 systems, both first 
principles calculation and angle-resolved photoemission spec- 
troscopy (ARPES) have demonstrated that a hole-pocket lies 
slightly below the Fermi level (Ep) at ambient P. With band 
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renormalization this hole pocket appears on the FS at some 
critical P, and the overall FS topology changes from only 
electron-pockets at P = to the coexisting electron and hole 
pocket as in pnictides. This topological FS transition induces 
a crossover from FS nesting along diagonal Qi — » (tt, 7r)] 
direction to along zone direction Q2 (tt, 0)/(0, it) which 
makes a pairing symmetry transition from nodeless d-wave to 
s ± -pairing. 

We explore these postulates via first-principle band struc- 
ture and pairing eigenvalue calculations within random-phase 
approximation (RPA). The first-principles calculation is per- 
formed within the atomic sphere approximation by using the 
TB-LMTO code. [021 [17) For a given P, the uniform volume 
contraction is evaluated using Birch-Murnaghan equation of 
state fT8l formalism (given below). The key here is to opti- 
mize the Se atomic position with respect to the Fe plane by 
minimizing total energy at each P. We find zs e — 0.3452c at 
P = GPa and z Se = 0.351c at P = 12 GPa. The obtained 
band structure in Fig.[TJc) shows that at P =12 GPa, the bands 
are renormalized by Z ~ 0.9 (see inset) which allows the 
hole-pockets to appear on the FS. Therefore the dominant FS 
instability arises along Q2 at a critical value of interaction U, 
and thus s ± pairing with large coupling constant commences. 
We also semi-quantitatively reproduce the phase diagram of 
the d and s ± -pairings within the conditions for best nestings 
along each channels. 

Formalism:- To grasp further insight into how P modifies 
electronic structure and also to enable adding correlations for 
pairing symmetry calculations, we derive a low-energy five- 
bands tight-binding (TB) formafism |[T2t to include pressure 
effect. The TB hopping integrals ti in a given crystal is de- 
fined as t = (Sk\V c \S&), V c is Coulomb interaction between 
lattice and electron, implying that ^(P)/^(0) ex V(P)/V , 
where U(P) and V(P) are the TB parameters and lattice 
volume at any given P P, and ^(0) and Vb are their cor- 
responding values at ambient P. For simplicity, we take 
U(P)/U(0) ~ V(P)/V « ra/ra* = Z, for all bands. From 
the value of Z or more strictly from the volume ratio, we 
can obtain the value of P using Birch-Murnaghan equation 
of statedl 
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where x = Vo/V = 1/Z, and P, B' are the bulk modu- 
lus and its P derivative. To match the optimal doping of 
the second T c dome with the experimental value we find 
B w B'P =67 GPa, which is close to the experimental 
value of B=64 GPa. [030 Since Z renormalizes the TB hop- 
ping parameters only, not the entire bands, Fermi level ad- 
justment is necessary at each value of P to keep the num- 
ber of electron constant. We compute the Fermi level self- 
consistently, beyond a simple rigid band shift approxima- 
tion, by integrating the density of state in the filled state as 
n = 1/N i=i-5 I du5{(jj—ZE l k ), where E l k is the disper- 
sion for i th band. Our finding of crucial importance is that P 
can enable significant FS topological change even with fixed 




FIG. 2. (Color online) (al-a3) TB bands and real-part of susceptibil- 
ity at zero energy at ambient P. (bl-b3) Same but at P = 12 GPa 
and Z = 0.9. (al), (bl) Self-consistently evaluated bands along 
high- symmetry lines in 1 Fe unit cell notation 1 12 1 at two P values. 
(a2), (b2) Corresponding FS topologies. (a3), (b3) Static susceptibil- 
ities plotted in the q = 0-plane. Arrows depict the leading nesting 
directions. 



number of electrons. In what follows, the Luttinger volume, 
Vl = l/NJ2k F i d(Ek F ) °f a U changes as a function 
of P but the ratio n = Vl/Z remains constant. This is due 
to the fact that, unlike in a Fermi liquid picture, here Z only 
renormalizes the dispersions, not the spectral function. 

Figure [2] reveals the evolution of electronic states at ambi- 
ent P and at P=12 GPa. At P = 0, the hole pockets lies 
at 60 meV below the Fermi level (Ef) at T -point, as is also 
seen in angle-resolved photoemission spectroscopy (ARPES) 
data. l20l In this case, the presence of electron pockets at X- 
points (in 1 Fe per unit cell notation) is well established.! 20- 
l22l With the uniform band renormalization by Z = 0.9 at 
P = 12GPa, according to Eq. 1, we find that two concen- 
tric hole pockets are fully formed on the FS as seen from 
Figs.|2[bl)and|2fb2). 

To determine the dominant nesting 'hot-spot', we perform 
multiband suscetibility calculation using full orbital overlap 
matrix-element effect. The many-body correction is incorpo- 
rated within RPA via explicitly including intraorbital interac- 
tion U, interorbital interaction V, Hund's coupling J > 
and the pair hopping energy J'. The full orbital compo- 
nents of interaction tensors U s for spin RPA susceptibility 
and U c charge RPA susceptibulity are defined in Ref. l23l 
and the details of the calculation can, for example, be found 
in Refs. El El El [25). We plot the static spin RPA sus- 
ceptibility (trace of the susceptibility tensor) at the two rep- 
resentative Ps in Fig. [2j a 3) for P = and in Fig. [2fb3) for 
P = 12 GPa. As expected, the dominant nesting in the for- 
mer case is aligned along the inter-electron-pocket direction 
[consistent with earlier calculation in Refs. lHTifT3ll l. On the 
other hand at P = 12GPa, the dominant nesting changes to 
Q2 — > (tt, 0) between electron and hole pockets, as obtained 
for many iron-pnictide superconductors . |[T0l l24l 

Next we evaluate the effective pairing vertex in the singlet 
channel for scattering between two FSs z, j within spin and 
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charge fluctuation exchange approximation! 25 ]: 
1 



P = 0GPa 



P = 12 GPa 
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The pairing strength for gap functions g(k) is then determined 
by solving the eigenvalue problem lP24l : 
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Here the line integrals over Q are performed on each FS 
loops, and vp is the Fermi velocity. The key point of Eq.[3]is 
that if the gap function g(k) possesses opposite sign at k and 
k\ mediated by a large peak in T at the 'hot-spot' q — k — k'. 
Based on this framework, we now study the evolution of A for 
g(k) = cos k x — cos k y for d x 2_ y 2 and g(k) = cos k x cos k y 
for s ± -pairing channels (in the 1 Fe per unit cell notation) as a 
function of P in Figs.[3]and[4](in 2 Fe per unit cell, the form of 
the pairing structure, i.e. g{k) transforms by the same unitary 
transformation as the FS so that the macroscopic properties 
such as nodeless and isotropic gap structure remains same in 
any notation! 13 , 26 ]). To understand the origin of the pairing 
symmetry transition in details, we study each intra- and inter- 
band components of A^ in Fig. [3] at two representative Ps. 
These results confirm our initial assumption. 

For ^2 -wave pairing at P = in Fig. |3ja) and at 
P = 12 GPa in Fig. (3fb), we obtain A > and < 0, re- 
spectively. This can be understood from the corresponding 
FS topological changes. At P = 0, the two electron pockets 
at (7r, 0) and (0, tt) posses opposite sign of SC gap, and this 
phase is supported by nesting along Qi — >> (tt, tt), as deduced 
in Fig.[^a3). l[TTUT3l So we get A 7 _ 7 > 0. But as the domi- 
nant 'hot-spot' changes to Q2 —> (0, tt) at P = 12 GPa, pos- 
itive A 7 _ 7 component is overturned by large negative A a _ 7 
and \(3- 1 components being supported by Q2 'hot-spot' with- 
out sign-reversal of g(k). Therefore, the total A becomes 
negative, making d x i_ y i pairing unstable at this P. For 
s ± — pairing, the situation is reversed in that A 7 _ 7 < at all 
values of P, but A a /£_ 7 > 0. It is important to note that due 
to the presence of two hole pockets a and /?, the total value 
of A becomes positive and large at some critical value of P, 
otherwise, s ± -pairing would have been favorable. 

In Fig. [4] we show the full P and interaction U depen- 
dence of the total pairing strengh A, and the corresponding 
calculated T c . The d-wave pairing survives upto P ~ 5 GPa, 
slightly less than the termination point of the first SC shown 
in Fig. [TJa). Of course, these quantitative consistencies rely 
strongly on the exact shape of FS topology and the value of 
bulk modulus B used in Eq. 1 . We focus on fitting B so that 
the optimum T c for the second dome match with the experi- 
mental value of P. With a separation of about 1-2 GPa, we 
find that the s ± — pairing channel appears abruptly for a large 
range of U. Although we obtain an optimum A as a function 
of P at which the FS nesting between a//3 to 7 is strong for 
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FIG. 3. (Color online) Evolution of different intra- and inter-band 
pairing eigenvalues, Xij , at two representative Ps (different columns) 
and pairing symmetries (different rows). Color coding on the FS in 
(al), (bl), (cl), and (dl) gives the variation and sign of SC gaps on 
each FS pockets. At P = (left column), the pairing eigenvalue 
arising from nesting between two electron-pockets is positive for d- 
wave, and negative for s^pairing, indicating that the former pairing 
is stable here. At P = 12 GPa, the strong negative value of 
and A/3- 7 , governed by same sign of SC gap for d-wave pairing con- 
nected by dominant 'hot-spot', yields total A < 0, and thus unstable 
pairing. For s ± -pairing at this P, despite the presence of several 
negative values of A, the dominant 'hot-spot' connects sign-reversal 
SC gap between two electron and two hole pockets, leading to total 
A > 0. If one of the hole-pocket disappears, s ± -pairing will become 
unlikely in this case. 



all values of U considered, it survives to a larger P range than 
the experimental data of K0.8Fe1.7Se2. However, for other 
samples within the same family, new data shows that T c in the 
second dome is very much P independent and survives up to 
P as large as 40 GPa measured so far. [27 ] 

We compute T c using standard McMillan's formula] 28 ] 
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A-/i*(l + 0.62A) 
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With a value of pseudopotential /i* = 0.1 and Debye tempera- 
ture uj d =U5 K, we obtain optimum T c -38 K and T c - 47 K 
for the first and second dome, respectively, for U = 1.2, 
which are close to the experimental values of 37 K and 48 K. 
It is obvious that the absolute value of T c depends on the in- 
teraction parameter, however, the ratio between the optimum 
value at two domes is always maintained. Interestingly, we 
find that the computed T c for the second dome is even flatter 
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FIG. 4. (Color online) (a) Evolution of total pairing eigenvalues for 
d-wave (solid circles) and s ± -pairing as a function of P for various 
values of U. The horizontal axis at the top of the plot gives the 
computed value of Z from Eq. 1 . The color shadings separates the 
negative and positing eigenvalue regions, (b) Computed values of T c , 
using Macmillan's formula in Eq. 4, are plotted as a function of P. 
The persistent of T c in the second SC region upto high P is observed 
in a number of materials other than K0.8Fe1.7Se2. (27) 



than that for A and thus agrees well with the new data. (27) 

Relating the glue function to the spin-resonance mode that 
appears in the SC state, we expect a resonance around uod ~ 
15 meV (~ 174 K) at the optimum P. The phenomena of 
spin-resonance in the SC state is well known lf8UTTl 24] and is 
essentially similar to the condition for obtaining positive pair- 
ing eigenvalue A as discussed above: Given that the sign re- 
versal of SC gap is connected by momentum transfer q a spin- 
resonance appears at an energy uj(q) ~ | A(fe) | + | A(fc + q) \ 
(RPA correction shifts the mode to a slightly lower energy). 
We compute the spin-resonance spectra at the same represen- 
tative P values as in Figs.[|and[3]for d-wave and s ± -pairing 
and the results are shown in Fig. [5] We immediately see 
that for d-wave pairing we get a nearly commensurate spin- 
resonance mode at uj/2A ~ 0.8, while it becomes incom- 
mensurate and shifts to higher energy for the s ± -pairing at 
higher P. Although, ARPES and INS measurements are dif- 
ficult to perform under P, neutron diffraction experiment can 
be done here to test the changes of Q vector as a function of 
P to verify our results. 

In conclusion, we present an analysis of pressure evolution 
of pairing interaction in A y ¥Q2- x ^2 family of superconduc- 
tors. We argue that observed two SC domes can be naturally 
explained by changes in FS topology driven by mass enhance- 
ments. At lower P, we expect no T -point hole pocket on the 
FS, and therefore <i-wave SC state. (29) Upon increase in P, 
we expect two hole pockets to develops at T-point, leading 
to pairing symmetry transition to s ± -pairing. Therefore, the 
natural outcome of our analysis is that there are two distinct 
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FIG. 5. (Color online) (a) Spin excitation spectrum along zone di- 
agonal direction as a function of energy for P = 0, and d-wave 
pairing, (b) Same but along zone boundary direction at P = 12 GPa 
and s ± -pairing, (c)-(d) Constant energy cuts at their corresponding 
resonance energy values. 



SC states, nodeless d-wave at smaller P and s ± -pairing at 
higher P that can be realized in these materials. To test these 
predictions, we suggest to use magnetic field dependent tun- 
neling spectroscopies to investigate gap changes. (30) Neutron 
scattering and neutron diffraction measurements would reveal 
the distinct pattern of spin resonance and nesting properties in 
these phases. 
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